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Abstract 

We generalize the Foppl-von Karman equations to an initially precurved sheet and present the 
underlying derivation. A geometrically computed moment of strain replaces the notion of bending 
moment and results in a geometric formulation of the theory of shells. As the curvature approaches 
zero, i.e., the sheet becomes flat, the new equations reduce to the classic Foppl-von Karman ones. 
The present theory solves the long-standing problem of formulating these equations for an a priori 
curved shell and applies, for instance, both to shell theory and to strongly curved biomembranes 
of cells as closed surfaces, exhibiting crumpling as the membrane thickness goes to zero. 
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Crumpling of thin membranes or sheets and the existence of auxetic (negative Poisson 
ratio^nmtCTials^have^ cha llenged phy sical imagination, and explanation, for quite a while 

. Figure Q shows a striking example stemming from 

se phenomena 



the biophysics of membranes [1 
occurring in thin sheets are realizations of the Foppl-von Karman (FvK) equations 
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2l| . an important case of the equilibrium equations of elasticity theory, incorporating 



nonlinearity, and aiming at describing strongly curved sheets of finite thickness h > 0. The 
FvK equations only apply, however, to a thin sheet that is planar to start with and is then 
bent strongly. Their derivation could not handle a sheet being precurved, although this often 
occurs in practice - as in Fig. 1, for instance. The latter problem has been tantalizing the 
physics of membranes ever since Foppl ^ and von Karman proposed their equations 
a century ago. 

The FvK derivation includes the assumption that the sheet in question is thin, but 
it allows for the possibility that deviations from the initial shape may be large: it goes 
to quadratic order in the normal deformation parameter. Because the sheet is thin, by 
assumption, it is still consistent to assume a linear stress-strain relationship within the 
sheet. Several studies have shown j2l| that this description makes physical sense for large 
deformation, up to and including sharp creasing - a regime well beyond what Foppl and 
von Karman originally had in mind. 

There is yet another assumption in FvK, which greatly limits its applicability to phys- 
ical systems of interest, namely, that the initial shape of the sheet to be strongly bent or 
even crumpled is planar. How would deformation of an initially curved, anisotropic surface 
interact with the pre-existing anisotropy? This and other fundamental questions of strong 
bending and crumpling of surfaces cannot even be asked in the context of the original FvK 
theory. We present here the generalization necessary to handle such questions. 

If one looks to shell theory, a notoriously difficult part of elasticity theory that deals with 
thin elastic surfaces [2^ , one finds that existing formulations are inadequate for the purpose 
of handling a priori curved sheets. Shell theory was developed in the context of mechanical 
engineering with practical applications in mind. The vocabulary of shell theory is entirely 
physical. One would not know from reading shell theory that there is already a highly 
developed mathematical language for dealing with surfaces and their deformations, namely 
differential geometry. In this paper we show how to use methods of differential geometry to 
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FIG. 1: Scanning electron micrograph of a 'ghost', a red blood cell (erythrocyte) voided of 
its cytoplasm through a small hole created by osmotic rupture of the cell (perhaps the dark area 
in the membrane). It consists of a closed cell membrane (plus reticulum) of thickness h ~ 10 nm; 
the cell's diameter is about 10 /im. The ridges give the membrane a crumpled structure. As was 
suggested by Lobkovsky J|, one can understand this "chaotic" structure through a sy stem of two 
coupled nonlinear partial differential equations due to Foppl 3] and von Karman |lJ | (FvK). One 
of the two differential operators of the form A 2 has a small prefactor scaling as h 2 where h is the 
membrane thickness; cf. (|14[) and (|24jl below. In the limit h — > it cannot be neglected but gives 
rise to crumpling - analogous to -R^A in hydrodynamics where it generates turbulence as the 
Reynolds number R e — ► oo. The FvK equations have been derived as a local deviation from the 
flat state. The cell membrane, however, was never flat and clearly shows that in biological physics 
this assumption often has no justification. Picture courtesy of T.L. Steck. The original has been 
contrast enhanced. 



formulate, and physically interpret, the equations of crumpling for a general sheet. It will 
be evident that this generality is necessary in describing deformations in which the change 
in curvature is comparable to pre-existing curvature. By way of motivation we call attention 
to a typical, biological, example in Fig. ^ an initially curved surface that has subsequently 
been wrinkled. 

We consider the equilibrium state of an elastic body occupying a region D C M. 3 and 
subject to a conservative force density / derivable from a potential U. The body is described 
by an elastic energy density E e \ quadratic in strain u^. Thus its equilibrium is characterized 
by 

= 6 [ dV (E cl + U)= f dV (a ij 5 Uij - f5u t ) (1) 
Jd Jd 

where the stress cr y ' is symmetric and linear in strain, and Sui is an arbitrary small defor- 
mation, with Suij a function of 8ui. Now we ask what this means in case the body is a 
thin sheet, meaning that the region D is a smooth surface Ai thickened along the normal 
direction by a small amount h, with Ai as the midsection. (Generalizations to such cases as 
h variable rather than constant, elastic constants dependent on position, etc., are straight- 
forward, and will not be emphasized in what follows.) Let the Euclidean metric of three 
dimensional space be gij, and let the first, second, and third fundamental forms of Ai be 
9af3i h a p, k a p, where Greek indices take values 1, 2. In particular the first fundamental form 
g a p is just restricted to Ai. We assume Ai is oriented by choice of a unit normal n. 

As is well known [25], there is an orthogonal coordinate system (q x ,q 2 ) on Ai such that 
the level lines of the coordinates are tangent to the principal curvature directions at every 
point (umbilic points may be coordinate singularities). Furthermore, in a thickened three- 
dimensional neighborhood of Ai this coordinate system may be extended to an orthogonal 
system (g 1 ,^ 2 ,^), where the third coordinate z is along the normal to Ai through {q x ,q 2 ), 
positive in the direction of n. Methods for computing geometrical objects within this frame- 
work have been given in [2^ . In particular, in these coordinates the metric g^ in Dm takes 
the form 

g^ = diag (5-11(1 - K X zf 522(1 - k 2 z) 2 , 1) . (2) 

Here diag(pn, 522) is the first fundamental form of Ai and (k 1: k 2 ) are the principal cur- 
vatures of Ai, functions of {q 1 ^ 2 ) but not z. The minus signs define the sign conventions 
relating k 2 ) and n. For example, if n is the inner normal on the sphere, then the princi- 
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pal curvatures are positive. There is still the freedom to choose (q 1 , q 2 ) such that at a given 
point P G M. the first fundamental form obeys gu t i = (722,2 = (at P). Throughout what 
follows )Q denotes a partial derivative d/dx a with respect to the coordinate x a following the 
comma. The only non-vanishing Christoffel symbols at P are now 

rL = = - l -g ll g n ,2 (3) 

r 22 = \g 11 922,i (4) 

and the corresponding ones with 1 <-> 2. 

A central ingredient in this approach is the observation that the two dimensional strain 
tensor on M. is related to the metric tensor by 

where 8g a p is the change in the metric tensor under a deformation V = {Su 1 , 5u 2 , if}), a vector 
field. Here the notation for the components of V has been chosen to suggest that tangential 
displacements Su a are small, but that the normal displacement if) is not necessarily small, 
reflecting the anisotropy of e.g. a biomembrane. The change in the metric q^fi under V , in 



turn, can be computed as a Lie- Taylor series, using the Lie derivative £y 
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Sg a p = (£vg)(d a , dp) + -(£ v (£ v g))(d a , d p ) + ... (6) 

Here £yg is the Lie derivative of the three-dimensional metric, with the result restricted to 
M. by evaluating on tangent vectors d a to M.. The result depends only on the restriction 
of V to Ai, and if V is given only on M., V may be extended in any smooth way for the 
purpose of the computation, for example the components of V could be independent of z. 
We therefore interpret Su^ as the strain in Dm, averaged across the thickness h. The result 
is, keeping terms linear in the deformation 5u a , and going to quadratic terms in if>, 

= (VJup + V (3 5u a + V a il)V(3ilj + il) 2 k a f 3 )/2 

-iph a p (7) 

where V Q is the covariant derivative on M.. (First derivatives of the function if> are, of course, 
just ordinary partial derivatives.) We illustrate the method by computing one component 
of the first Lie derivative. The vector field V is the first order differential operator 

V = 5u 1 d l + 5u 2 d 2 + ipd z . (8) 



Then we have, evaluating at P for simplicity, and noting that z = there, 

{£ v g)(d 1 ,d 1 ) = Vg(d 1 ,d 1 ) + 2g([V,d 1 ),d 1 ) 

= Su 2 g ni2 - 2ipn 1 g 11 + 2g n 5u\ 
= 2(Su 1A + <Ju 2 rii) - 2if}g n Ki 

and similarly for other components. Since the result is a tensor and coincides with V Q <5«/3 + 
V p8u a — 2iph a p in this coordinate system, it must be this tensor. The second derivative 
computation is similar. 

Since the three-dimensional body D^i is thin, we assume the strain is at most linear 
in z. But when Dm is bent, it is relatively compressed on one side and extended on the 
other, so it is clear that there is a physically important strain that is linear in z. From the 
geometrical meaning of the second fundamental form h a p as the rate of rotation of the unit 
normal n, it is clear that this contribution is —z5h al 3, where 8h a p is the change in the second 
fundamental form of M. under the deformation. 

Unlike the first fundamental form, the second fundamental form is not the restriction 
of a tensor field in three dimensions to Ai, and so its change under deformation must be 
computed in a different way. For example, the second fundamental form of the varied surface 
Ai' is, up to a factor —2, the Lie derivative of g with respect to the unit normal on Ai', 
according to (J2J), and there are also more classical ways to compute it j^jj]. The result is 

$ U al3 = -zSha/3 = -z(S/ a V^ - ^k afi ) . (9) 

Here we have only gone to linear order in ip, as the bending term is usually small in any 
case, and the second order terms in ip are complicated. Again we illustrate the method. Let 
the equation of the varied surface Ai' be z = ip, the result of deforming A4 by the vector 
field U = ipd z . Tangent vectors to M! are 

X li = d li + [d li ,U]=d^ + ^ >li d z . (10) 

The unit normal n' on M! is orthogonal to these and normalized, hence to first order in if) 

n' = -g 11 ^ A d 1 -g 22 tP >2 d 2 + d z . (11) 
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Then computing one component, at a point over P for simplicity, we have 



(£ n , )(Xi,Xi) = n f g(X 1 ,X 1 ) + 2g([X 1 ,n'},X 1 ) 

= -^,2 5 ,22 5'll,2 + 25ii(l - £Kl)(-/Ci) 

We evaluate at z = if), finding 

{£n>g){X 1} Xi) = -2 ffll K! + 2^HK? - 2(^,11 + ^, 2 r?i) 

and similarly for other components, giving the result above. Going to higher order in if> 
would only require solving to higher order in (|11|) as was done, for example, in 2^|. 
As stress is linear in strain, it too will be linear in z: 

o" = at + zof . (12) 

Putting 5u af 3 = 5u£p + Su^g from (JJJ), 0, and a a/3 from (fT2j) into (JTJ), and integrating z 
from —h/2 to /i/2 we obtain 

0= / dAy/g[ha%6v% + ^ogfo%-h(J a 6u a + f'1>)] 
Jm i/ 

where g := det(g a g). Finally, integrating by parts and ignoring boundary terms (if A4 is 

a closed surface, for example), and recognizing that the variations are arbitrary, we obtain 

the equations of equilibrium 

= V^Jf + f a , (13) 

+ hatiKp + V Q V^ - k a pif)) +h(f z - V a if)f a ) , (14) 

a nonlinear system of equations for deformation (u a , if>) in response to the force field /. 

It is interesting to notice that one effect of keeping second-order terms in if) in (j2J) is to 
correct the second fundamental form h a p in the second parenthesis of (|14j) for the change 
in curvature due to the normal displacement if> - compare which contains this same 
expression, arrived at differently. Thus one effect of the nonlinearity is to replace the original 
h a p on M. by h a g(if>) in the balance of normal stress. This can be a large effect, because 
the tangential stress can be large, and it is frequently the o"^/i Q/ g term which is most 



important in balancing applied normal stress. If a* = £<? a ^, for example, corresponding to 
surface tension £ in Dm, then that term is the normal stress EH, where H = («i + K2) is 
the mean curvature. 

The insight of Foppl and von Karman can be appreciated by taking the special case of an 
initially flat A4, so that h a p = k a p = 0. Then, taking Cartesian coordinates, the covariant 
derivatives are ordinary derivatives, and the equilibrium equations (|T3|l and (fTlj) reduce to 

= dp** + r , (15) 

= ^{dadpaf) + ha*(d a d^) + hf z . (16) 

Equation (fTKJ) is essentially the second Foppl-von Karman equation. Its middle term, rep- 
resenting the normal stress due to a*, is absent without the nonlinear term in (|7Jl. leaving 
just bending stress to balance normal stress. The corresponding linear theory, which this 
theory was designed to correct, greatly underestimates the strength of the membrane to 
resist normal stress. As we see in (114)1 . the contribution of the nonlinear term persists in 
the case of a curved A4, but its effect is less dramatic. 

It is worth examining other aspects of the flat M. theory to see what else persists in the 
more general case. The usual linear relation between stress and strain, which we have not 
yet invoked, should still hold in general, 

E 



a 



a/3 



ig^du^ + ae^du^) (17) 



1 - a 2 

where E is Young's modulus, a is Poisson's ratio, and e af3 is the antisymmetric tensor with 
e 12 = 1/y/g and g = det(g a/ 3). In particular, 

E 



a 



* - (g a ^ u 5u* + ae^e^Su B ) (18) 



with 5u B given in Q. Equivalently, strain and stress are related by 



and in particular 



In the flat case 



5u a p = —{g^g^o^ - ae^e^a^) (19) 



5u Z = -^(gan9PvVZ - Ota^PvOZ) • ( 20 ) 



dfpt = (21) 
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implies is derivable 



29j from an Airy potential \i 



a 



ap 
M 



e a ^d,d uX (22) 



eliminating one variable. In this case, putting (}2"2*|) into the right side of ()2(Jj) and from 
(J7J) into the left side of (|2Uj). and applying the operator e aK ep\d K d\ to both sides, we find the 
first Foppl-von Karman equation 

-det(9 a ^) = iA 2 X (23) 

and putting (|2*2*|) and (fTHj) into (JTSJl we obtain the second Foppl-von Karman equation 

= -k c AV + h*?v<? v (dtd v x){d a dffll>) + hf z . (24) 

Here k c = Eh 3 /[12(1 — a 2 )] is the bending rigidity, with K c /h oc h 2 . Tangential displace- 
ments 5u a being small, det(d a dpif)) in (|2*3j) equals the Gaussian curvature K = kik 2 of the 
membrane surface M. to fair approximation. Bending in two orthogonal directions, which 
is what bending in general boils down to, implies j^O] IK ^ 0. 

The general equations (fT3*|) and (|14jl solve the long-standing problem of describing dy- 
namic equilibrium of a precurved sheet under strong deformations. They can be widely used 
in shell theory and, for instance, to analyze crumpling 0,0,0] of a naturally precurved cell 
membrane, such as the one in Fig. Unfortunately the above program for an originally flat 
surface cannot get started in the general, curved case. The non-commutativity of covariant 
derivatives means there is no Airy representation to eliminate (|13j) even if f a = 0. The equi- 
librium equations, using all the tricks known for less general situations, remain Eqs. (fT3*j) 
and (j!4|) . with o~ a " given by (|17|). Through the second and third fundamental form (h a p) 
and (k Q p) and the covariant derivatives they explicitly show how curvature must be taken 
into account. Evaluating the consequences of these equations, especially that of taking the 
limit h — > 0, will be a true challenge to membrane physics for some time to come. 

The authors thank Ted Steck for providing Fig. 1. J.L.v.H. gratefully acknowledges 
constructive discussions during a stimulating workshop at the Aspen Center for Physics, 
where it all started. M.A.P. thanks the SFB 563 for support while this work was done at 
the TU Munich. 
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